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ABSTRACT 

Using RXTE and ASCA data, we investigate the roles played by occultation and 
absorption in the X-ray spin pulse profile of the Intermediate Polar PQ Gem. From 
the X-ray light curves and phase-resolved spectroscopy, we find that the intensity 
variations are due to a combination of varying degrees of absorption and the accretion 
regions rotating behind the visible face of the white dwarf. These occultation and 
absorption effects are consistent with those expected from the accretion structures 
calculated from optical polarisation data. We can reproduce the changes in absorber 
covering fraction either from geometrical effects, or by considering that the material 
in the leading edge of the accretion curtain is more finely fragmented than in other 
parts of the curtain. We determine a white dwarf mass of ~1.2 using the RXTE data. 

Key words: binaries: close - stars: individual: PQ Gem - stars: magnetic fields - 
novae, cataclysmic variables - X-rays: stars - accretion 



1 INTRODUCTION 

Magnetic cataclysmic variables (MCVs) can be split into two 
groups: those in which the magnetic field of the accreting 
white dwarf is strong enough (<; 10MG) to synchronise its 
spin period with that of the binary orbital period - the po- 
lars, and those with a magnetic field insufficiently strong to 
achieve this synchronisation - the intermediate polars (IPs). 
The MCV PQ Gem is unusual in that it exhibits characteris- 
tics of both groups: it shows a strong soft X-ray component 
(kT ~50eV) (Duck et al. 1994), it is polarised in the op- 
tical/IR wave-bands (Potter et al. 1997; Piirola, Hakala & 
Coyne 1993) and has an estimated magnetic field strength 
of ~8-21 MG (Vath, Chanmugam & Frank 1996; Potter et 
al. 1997; Piirola et al. 1993) - all of which are character- 
istics of polars. On the other hand it shows a spin period 
of 833.4 sec (Mason 1997) and an orbital period of 5.19 hrs 
(Hellier, Ramseyer & Jablonski 1994) which are typical of 
IPs. PQ Gem can therefore be thought of as the first true 
"intermediate" polar (Rosen, Mittaz & Hakala 1993). 

PQ Gem has been observed using several X-ray satel- 
lites (ROSAT, Gmga, ASCA & RXTE). The X-ray light 
curves show a prominent modulation on the spin period, 
in particular, a pronounced dip in the light curve which is 
thought to be due to an accretion stream obscuring the main 
emission region on the surface of the white dwarf. Mason 
(1997) made a study of the then available X-ray data to de- 
termine an accurate ephemeris for PQ Gem based on timings 
of the dip. However, a detailed study of the spectral infor- 
mation contained in the ASCA data was not undertaken. 



This paper is primarily targeted at reaching a greater 
understanding of the interplay between the emission sites 
and absorption which produces the observed modulation of 
the X-ray light curves. This is achieved through analysis 
of the ASCA spectral data with supporting evidence from 
the hard X-ray RXTE light curves. The mass of the white 
dwarf is calculated from the RXTE spectral data using the 
stratified accretion column model of Cropper et al. (1999, 
subsequently CWRK). A fuller report is available in James 
(2001). 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 ASCA 

ASCA was launched in 1993 carrying two X-ray CCD cam- 
eras (SIS) and two imaging gas scintillation proportional 
counters (GIS), and operated until 2000 (Tanaka, Inoue & 
Holt 1994). The SIS detectors covered the energy range of 
0.4 keV to lOkeV with an energy resolution of 2% at 5.9 
keV. The GIS detectors had an energy range of 0.7 keV to 
10 keV with energy resolution of 8% at 5.9 keV. Above 8keV 
the GIS had a greater effective area than the SIS. 

Details of the observations of PQ Gem are given in Ta- 
ble The SIS detectors were configured to faint data mode 
2-CCD (high bit rate) and 1-CCD (medium bit rate) clock- 
ing modes and the GIS detectors were configured to PH- 
mode. It is not possible to merge data from the different 
types of detector or from different SIS clocking modes with- 
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Date 


Instrument 


Bit 

rate 


Integration 
time (s) 


Average 
cts/s 


1994/11/04 


SISO 


high 


26786 


0.35 






medium 


36910 


0.38 




SISl 


high 


24732 


0.29 






medium 


34734 


0.31 




GIS2 


high 


27771 


0.25 



Table 1. Details of the ASCA observations and data used in the 
analysis for this paper. 



Date Instrument Integration Average 

time (s) cts/s 

1997/01/27-30 PCA 51216 13.9 



Table 2. Details of the RXTE PCA observations and data used 
in the analysis for this paper. 



out loss of information. Hence the end product of the data 
selection and reduction was several data files with effective 
exposures ranging from 24.7-36.9ksec. 

The data selection and reduction were carried out fol- 
lowing recommended procedures (details in James 2001). 
The observation was made early in the ASCA mission prior 
to the onset of the degradation of the CCDs due to radiation 
damage (Yaqoob et al. 2000). Hence, the background spectra 
for SIS were created using the appropriate blank-sky event 
lists (available as part of ASCA calibration data) together 
with the data selection criteria corresponding to each event 
list. Background spectra for CIS were created in a similar 
fashion using blank-sky event lists. Background light curves 
were created from a source-free region of the same observa- 
tion and subjected to the same screening criteria as for the 
source light curves. 



3.1 Period Analysis 

The variability in the RXTE and ASCA light curves was 
analysed using a standard discrete Fourier transform (DFT) 
code (Deeming 1975, Kurtz 1985). The equation, zo = 
— ln[l — (1 — Po) 1 ], (Scargle 1982) was used to determine 
the power level, zq, above which an amplitude peak would 
be spurious for a small fraction, po, of the time, where N 
is the number of frequencies examined. A 90% confidence 
(po = 0.1) was used to determine the signal-to-noise level 
above which the amplitude of a power peak was judged to 
be significant. The noise level of the periodogram was taken 
to be the mean of the amplitude spectrum after it had been 
prewhitened with frequencies found to be significant (Table 
B|). The error in each such frequency was taken to be the 
standard deviation, a, from a least squares fitting. 

The ASCA data showed three significant periods, 
833.46s, 416.678s and 277.763s which are consistent with 
spin ephemeris of Mason (1997) and its first two harmon- 
ics. The RXTE data showed a significant period of 833.54s. 
The spin period (833 sec) is consistent with the spin down 
rate predicted by the ephemeris of Mason (1997). The com- 
plete results are given in Table ^ for those signals identi- 
fied as significant at the 90 percent confidence level. The 
amplitude spectra, and window functions for the RXTE 
and ASCA are shown in Figure [j], along with the spectra 
prewhitened with the frequencies of the significant power 
peaks and the expected positions of the orbital, beat, spin 
frequencies and spin harmonics. The residual power peaks in 
the prewhitened RXTE periodogram are the 3rd harmonic 
of the spin period (> 68% confidence level) and the 1st and 
2nd harmonics (<68% confidence level). 

Using an orbital period of 5.19hrs (Hellier et al. 1994) a 
beat period at 14.54 min may be expected. However, there 
is no evidence for a significant amplitude at the frequencies 
corresponding to the orbital or beat periods in either the 
RXTE or the ASCA data (Figure |). This indicates that the 
accreting material must go through a disc since all orbital 
information will then be lost. 



2.2 RXTE 



The Rossi X-ray Timing Explorer (RXTE) was launched on 
30th December, 1995. A review of the available instrumen- 
tation and the scientific results of the first years are given 
in Swank (1998). 

Details of the RXTE observations of PQ Gem used in 
this paper are given in Table [| Only the top Xenon layer 
was selected from the five PCA instruments, so as to im- 
prove the signal-to-noise ratio. The standard RXTE proce- 
dures were used to select data (details in James 2001). The 
background was modelled using the latest faint background 
model appropriate for the epoch of the observation, i.e. files 
faintl7_e03v03.mdl and faint240_e03v03.mdl. 



3.2 ASCA light curves 

Light curves of these data were first presented in Mason 
(1997) and also in James et al. (1998). They are presented 
here for completeness and comparison with those from 
RXTE. Light curves using data from the SISO detector for 
energy bands 0.7-1. OkeV, 1.0-2. OkeV and 2.0-4.0keV, 4.0- 
lO.OkeV and from the total bandpass (0.7-10.0keV) folded 
on spin period are shown on the left-hand side of Figure ^. 
The light curves were not folded on either the orbital or beat 
periods due to the lack of amplitude peaks found at thes e 
frequencies in the power spectrum of the ASCA data (§3.1). 



3 LIGHT CURVES 

The quadratic spin ephemeris of Mason (1997) was used 
throughout the analysis of the light curves which were he- 
liocentrically corrected. 



3.3 RXTE light curves 

Light curves were extracted for the energy bands 2. 0-4. OkeV, 
4.0-10.0keV, 8.0-25. OkeV and from the total bandpass, 2.0- 
25keV, which were then folded on the spin period (right- 
hand side, Figure [|). The first two energy bands can be 
directly compared with similar plots using the ASCA data. 
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Figure 1. Periodograms for the ASCA 0.5-10.0keV light curve (upper half) and RXTE 2.0-25.0keV light curve (lower half) of PQ 
Gem. In each half the top plot shows the window function, in the middle plot, the amplitude spectrum, and the bottom plot shows the 
spectrum prewhitened with the frequencies of the significant peaks. The vertical bars (from left to right) mark the positi on o f the orbital, 
beat and spin (w) frequencies plus the 1st, 2nd and 3rd harmonics of u), respectively. The relevant details are given in §|3.l| and Table H. 
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Figure 2. Background-subtracted light curves of PQ Gem using ASCA SISO and RXTE PCA data folded 
on the spin period. On the left hand side are the ASCA light curves for energy bandwidths 4.0— lOkeV (top), 
2.0-4.0keV (second from top), 1.0-2.0 keV (middle) and 0.7-1. OkeV (second from bottom) 0.5-10.0keV 
(bottom). On the right hand side are the RXTE light curves for energy bandwidths 8. 0-25. OkeV (top), 
4.0-10.0keV (second to top), 2.0-4.0keV (third to top), 2.0-25keV (bottom). 



3.4 Comparison between the ASCA and RXTE 
light curves 

The RXTE PCA instrument overlaps the ASCA SIS instru- 
ments across the 2-10keV energy band. Hence this energy 
band from both instruments was used to compare the light 
curve modulation of the two observations (Figure []). The 
light curves in these 2 bandpasses were found to be very sim- 



ilar, each with two maxima separated by a minimum at spin 
phase <f) ~ 0.0 (referred to as "the dip", Mason et al. 1992) 
and a second minimum at (f> ~ 0.6. In the 4. 0-10. OkeV light 
curves the minima are relatively shallow in both the ASCA 
and the RXTE light curves, but the 2.0-4.0keV light curves 
show well defined "dips" at (f> — 0.0, with a broader mini- 
mum at (f> ~ 0.6. There is some evidence that the amplitude 
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Satellite 


Instrument 


Passband 
keV 


Frequency 
mHz 


Period 

s 


S/N 


ASCA 


SISl 


0.7-10.0 


1.19981(25) 


833.46(17) 


8.0 








2.39994(32) 


416.678(56) 


7.1 








3.60019(32) 


277.763(25) 


6.6 


RXTE 


PCA 


2.0-25.0 


1.199703(81) 


833.540(56) 


12.8 



Table 3. The details of the significant power peaks found from period analysis of ASCA SISl data and RXTE PCA data. The signal- 
to-noise ratio was taken as the amplitude of power peak to the mean of the prewhitened spectrum. 



of the variation in the RXTE folded light curve is less than 
that in the ASCA SIS data. The shape of "the dip" in the 
2.0-4.0keV band is "V"-shaped in the RXTE compared to a 
more "U" -shaped in the ASCA light curve and its depth ap- 
pears greater in the latter. This may be due to a higher level 
of absorption at the earlier epoch of the ASCA observation. 



3.5 RXTE hardness ratios 

Plots of the 8.0-25.0keV/2.0-8.0keV and 8.0-25. OkeV/2.0- 
4.0keV hardness over the spin period are shown in Figure ^| 
The lower limit of the 8.0-25.0keV energy band was chosen 
so as to be well above the main absorption edges. A higher 
value of the ratio indicates those positions in the spin cycle 
which are more affected by absorption. These are consistent 
with the dip at <j> = 0.0 being due to absorption. 



4 SPECTRAL ANALYSIS 

Because the energy resolution of the RXTE PCA (18% at 
6.0keV) is much poorer than that of the ASCA SIS (2%), 
and does not extend to lower energies, the spectral analy- 
sis was restricted to the ASCA data. During the spectral 
analysis the medium and high bit rate SIS data from both 
instruments and the GIS data were linked through the mod- 
els and fitted simultaneously. 



4.1 Integrated Spectrum 

The integrated spectrum was modelled using the MEKAL 
code for an optically thin emission model (Mewe, Kaastra & 
Liedahl 1995) for the hard X-ray spectrum and a blackbody 
model for the soft X-ray continuum. The temperature of the 
former was difficult to constrain, a problem which was sim- 
ilarly experienced by Duck et al. (1994) using Ginga data. 
Therefore the temperature was fixed at 20keV in line with 
the results from their analysis (we tried fixing the tempera- 
ture at other values, but this did not have a significant effect 
on the results). The temperature of the blackbody compo- 
nent was not well constrained and therefore this spectral 
component was fixed at 55eV in line with that found with 
ROSAT data (Duck et al. 1994) due to its greater sensitivity 
to this spectral range. The fluorescent iron line at 6.4keV was 
modelled using a Gaussian component, the width of which 
was fixed at 0.05keV. A single homogeneous photo-electric 
absorber gave a poor fit (xt = 2.0), therefore an inhomoge- 
neous (partial covering) photoelectric absorber was added to 
the model. A good fit to the data was achieved in this way 




RXTE PCA 



8.0 - 25.0keV / 2.0 - 4.0keV 



RXTE PCA 



8.0 - 25.0keV / 2.0 - 8.0keV 



1.0 
Spin Phase 



Figure 3. Plots of the 8.0-25. 0keV/2.0-8.0keV (below) and 
8.0-25.0keV/2.0-4.0keV (above) hardness ratios folded on the 
spin period. The errors are la. 



(xt = 1-07). The values of those parameters allowed to vary 
during the fitting are given in Table ^ and the integrated 
spectra are plotted in Figure El 



4.2 Phase-resolved Spectroscopy 

Even though a good fit to the integrated spectrum was ob- 
tained, the results from the 8.0-25.0keV/2.0-8.0keV and 
8.0-25. 0keV/2.0-4.0keV hardness ratios using the RXTE 
PCA data indicate the possibility of variation in the ab- 
sorption during the spin cycle (Figure H, §3.5). These two 



results imply that variability in both the absorption and the 
normalisation of the hard X-rays over the spin cycle could 
be the cause of the observed modulation in the spin-folded 
light curves. In order to investigate this further we divided 
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Nff(l) N H (2) (fraction) BB MEKAL Gaussian Flux(observed) \l (dof) 

10 22 cm -2 10 22 cm~ 2 norm: norm: norm: 2-10keV 

10" 3 10~ 2 1CT 5 10 _11 erg cm" 2 s" 1 

0.12±°;° 2 8.3+?j (0.53+°;°^) 2.0+ 2 ' 3 1.7lt°;^ 2.12 1.07 (1070) 



Table 4. Results from spectral analysis of the ASCA integrated spectrum of PQ Gem including the error range at 90% confidence 
level. The optically thin plasma and blackbody temperatures were fixed at 20keV and 55eV respectively. Njj(l) corresponds to the 
homogeneous absorber whereas N/j/(2) applies to the inhomogeneous (partial covering) absorber; the flux is not corrected for absorption. 
The model was fitted simultaneously to the medium and high bit rate data from the SIS1 and SIS0 instruments and the high bit rate 
data from the GIS2 instrument. 




12 5 10 

channel enengy (keV) 



Figure 4. Integrated spectra of PQ Gem. Data from the 
medium and high bit rate data from the SIS1 and SIS0 instru- 
ments and the high bit rate data from the GIS2 instrument 
with the fitted model (Table 0). 



the spin period into 10 equal time bins and extracted spectra 
for each using data from both the SIS and GIS2 instruments. 
During all subsequent spectral fitting the medium and high 
bit rate data from both SIS and the GIS2 instruments were 
linked through the models and fitted simultaneously. 

The reference model was derived from 0.4 < <j> < 0.5, 
since this phase had the greatest number of counts. The 
model resembled that obtained from the integrated spec- 
trum with the exception of the blackbody component whose 
temperature was significantly higher (best fit=83eV). In this 
more detailed analysis it was considered to be more appro- 
priate to fix the temperature of the blackbody component 
of the other phase resolved spectra at this value. The best 
fit model comprised of three emission components: an opti- 
cally thin plasma emission component with a temperature of 
20keV (fixed), a blackbody component with a temperature 
of 83eV and a Gaussian at 6.4keV to model the fluorescent 
Fe emission line (fixed width of 0.05keV). Two photoelectric 
absorption components were required to model the atten- 
uation of the X-rays: an inhomogeneous (partial covering) 
neutral absorber and a homogeneous neutral absorber. The 
latter had a best fit of 6.8 x 10 20 cm -2 : this component was 
fixed at this value in the other phase intervals. The metal 
abundance had a best fit of 0.34 solar - again this was fixed 
at other phase intervals. We now allowed the normalisation 
of the X-ray emission and the partial covering absorption 
component to vary during the remaining 9 phase bins of 



the ASCA data. We also allowed the normalisation of the 
blackbody and the Fe fluorescent emission line to vary. 

The resulting best fit parameter values from this analy- 
sis are shown in Table |^. The variations in the normalisation 
of the optically thin plasma, the column density and partial 
covering fraction are plotted in Figure |f| together with the 
folded light curve for the 0.5-10.0keV energy band of the 
ASCA data. It can be seen that the column density does 
not vary significantly over the whole spin cycle. The nor- 
malisation parameter of the hard X-ray component shows a 
maximum over 0.2 < <f> < 0.4, decreasing to a minimum at 
(f> — 0.7. The covering fraction of the partial absorber com- 
ponent shows a steady decline from its maximum value at 
= 0.0 to a minimum value at (j> = 0.5. This pattern and 
phasing closely follows that of the RXTE 8.0-25. OkeV/2.0- 
8.0keV hardness ratio, shown in Figure 3. 

Relating this back to the main features of the modula- 
tion of the light curve shows that the maximum countrate 
(0.1 < (f> < 0.4) is due to a combination of increasing nor- 
malisation coupled with a decreasing level of absorption. It 
would appear that an increase in the normalisation factor 
is probably the main cause of the second maximum (0.75 
< 4> < 0.85) of the light curve, but a contribution from a 
decrease in the covering fraction may also play a part. The 
secondary minimum (0.55 < (f> < 0.75) is mainly due to 
a lower normalisation factor. Finally, "the dip" can be ac- 
counted for by the maximum value of the covering fraction, 
although, again, a contribution from a slight decrease in the 
normalisation factor cannot be ruled out. 



5 MASS OF THE WHITE DWARF 
5.1 Procedure 

In recent years much work has been done in determining the 
shock temperature in magnetic CVs (e.g. Cropper, Ramsay 

6 Wu 1998, Fuijmoto & Ishia 1996). One of the main goals 
of this work has been to estimate the mass of the white dwarf 
since these two parameters are closely linked. To estimate 
the mass of the white dwarf in PQ Gem we use the emission 
model of Cropper et al. (1999) (CWRK) to fit data from 
the RXTE PCA detector (since it has an extended high 
energy range - the temperature of the hard X-ray component 
was poorly constrained using the ASCA data because of the 
lOkeV upper limit). This model was developed to fit high 
quality X-ray spectra appropriate to the conditions found 
in the accretion sites of MCVs. It is a more complex model 
than the generalised optical thin plasma code we used for 
the phase-resolved spectroscopy. 
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Spin Phase N ff (2) Covering BB MEKAL Gaussian: Flux(observed) xt ( d °f) 

10 22 cm -2 Fraction Norm: Norm: Norm: 2— lOkeV 

lCT 4 10~ 2 10~ 5 10 -11 erg 

— 9 — 1 

cm s 



0.0-0.1 


7 


4+i.s 

— 1.5 





on+0-02 
—0.02 





48+V^ 
"—0.48 


1 


88+° 

u — 


21 
18 


7.96+? 

' *^ — 6 


94 
33 


1.86 


1 


09 


(149) 


0.1-0.2 


8 


" — 2.2 





70+°°:? 

— 0.04 


1 


62+°" 
u — 1.60 


2 


10+° 

— 


2(i 
22 


11. 9+? 

^—5 


20 
64 


2.17 


1 


00 


(198) 


0.2-0.3 


5 


— 1.3 





kq+0.03 
—0.03 


1 


24+0.51 
—0.54 


2 


22+?, 
"— 


20 
16 


12. 7+? 


20 
6!) 


2.59 


1 


11 


(248) 


0.3-0.4 


5 


7 +2.2 

-1.5 





kc+0.03 
—0.04 


1 


63+°?? 

—0.64 


2 


21+9 
— 


21 
17 


7.77+1 


53 
21 


2.55 





91 


(245) 


0.4-0.5 


8 


-1+3.9 
-2.7 


o 


Aj+O-OG 
^'-0.05 


1 


ic+0.55 
±u -0.38 


2 


18 +0 


28 
2(1 


8 00+ 3 
o.uu_ 5 


50 
Oil 


2.49 


1 


12 


(253) 


0.5-0.6 


3 


O+2.0 
8 -1.5 





43+0.04 
^°+0.04 


1 


A7 +0.55 
*'-0.43 


1 


791° 


17 
14 


7.02+. 2 


68 
94 


2.23 


1 


08 


(241) 


0.6-0.7 


6 


q+4.1 

a -2.6 





54+0.05 

d ^-0.06 





rro+0.48 
OO -0.35 


1 


691° 


23 
17 


11.7+J 


30 
42 


1.96 





98 


(202) 


0.7-0.8 


1 


+1.9 
u -1.3 





c- 9 +0.03 
J -0.04 





ck+0.49 
DO -0.36 


1 


61±S 


13 
12 


6.69+J 


51 

8!) 


1.96 





99 


(209) 


0.8-0.9 


8 


i+3.9 
-2.7 





c n +0.05 
OU -0.05 





C9+0.47 
d -0.46 


1 


991° 


20 
18 


4.21±4 


98 
(Hi 


2.23 


1 


04 


(231) 


0.9-0.0 


4 


0+1.5 
8 -l.l 





co+0.03 
°°-0.03 


1 


49+0.48 
* z -0.55 


1 


871° 


l(i 
12 


11.8^ 


20 
24 


2.2 





95 


(218) 



Table 5. The variable parameter values from spin phase-resolved spectral analysis over 10 phase-bins. The reference model (0.4 < <f> < 
0.5 ) was simultaneously fitted to the medium and high bit rate data from both of the SIS instruments and high bit rate data from the 
GIS2 instrument on the ASCA satellite. The normalisation and absorption parameters only were varied from those of the reference model. 
The blackbody, bremssstralung normalisations and absorption include the error range at the 90% confidence level. The blackbody, Mekal 
and Gaussian normalisations are standard for the relevant xspec models. Also included are the x 2 and observed flux (i.e. not corrected 
for absorption) for the 10 phases. 



The parameter, e s (the ratio of the cyclotron to 
bremsstralung cooling at the shock), was fixed at 0.5 which, 
from figure lb of Wu, Chanmugam & Shaviv (1995) corre- 
sponds to a magnetic field strength of ~ 15MG (Piirola et 
al. 1993, Vath et al. 1996, Potter et al. 1997). We stratified 
the volume between the surface of the white dwarf and the 
shock front (Cropper et al. 1998) into 100 levels. The free 
parameters during the spectral fitting were the specific ac- 
cretion rate in g cm _2 s _1 , and the mass of the white dwarf. 
A viewing angle is also required to scale the reflection com- 
ponent in the model. The abundance was allowed to vary 
from a solar profile in the final stages of the fitting. 



5.2 The Results 

The RXTE spectral fitting was restricted to the spin phase 
range of the data to 0.1 < (f> < 0.5 which corresponds to that 
portion of the spin cycle where a single emission site is visible 
(Potter et al. 1997). This phase range also excludes the dip 
centered on 4> =0.0 which from ROSAT and ASCA data, 
has been shown to be due to absorption. Models including a 
homogeneous absorber components together with a CWRK 
emission component were found to be sufficient to produce 
an acceptable fit (xl ~0.89) to this data using an estimated 
average viewing angle of ~60°. This gave a white dwar f mass 
of 1.21M (Table |). We discuss this further in ^6. 



6 DISCUSSION 



In section 4.2 we find that the spin phase modulation can be 
modelled well by a variation in the partial covering fraction 
of a neutral absorber at low energies and in the normalisa- 
tion at higher energies. Before we discuss this further, we 
briefly examine whether the modulation at both high and 
low energies can both be explained by a tall accretion shock. 

In the case where the shock has significant vertical ex- 
tent, the lower (cooler) part of the shock could be obscured 
by the limb of the white dwarf at certain spin phases. The 
light curves will then have a larger amplitude at low ener- 
gies than at high energies (e.g. Allan, Hellier & Beardmore 
1998). The height of the accretion column, H, can be esti- 
mated from the relationship 



H 



5.45 x 10 8 M 1 - 6 1 /-2M^ J R^ 2 D 



(Frank, King & Raine 1991). Using the lower limit of the 
fractional area and the corresponding M derived in § |6.4| (cf 
Table (?]) gives an upper bound of H=5Ax 10 7 cm or 0.14 
Kwd- However, the observational evidence suggests that PQ 
Gem has a significant magnetic field (§jjj) which implies that 
the shock height will be lower than this (e.g. Cropper et 
al. 1999) due to cyclotron cooling. This suggests that this 
mechanism is not the cause of the variation. However, the 
accretion regions are likely to be sufficiently structured (e.g. 
Potter et al 1997) so that such a scenario cannot be ex- 
cluded. On the other hand, because the spectral variations 
can be well modelled by a variation in covering fraction and 
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Satellite detector angle N^f co iij(l) fa re WD mass x 2 (dof) 

° 10 22 cm -2 gcm" 2 s" 1 10 7 cm M & 

RXTE PCA ~60 5.2 2.3 9.5 1.21 (1.16-1.28) 0.89 (46) 



Table 6. Results from fitting the CWRK model to the RXTE PCA spectrum (0.1 < <j> < 0.5). Details of the spectral fitting process are 
given in §H 



WD mass \t (dof) unabsorbed flux luminosity M fractional 

(0.001-100.0keV) area 



(g cm 2 s 1 ) 






(ergs cm 2 s x ) 


(ergs s 1 ) 


(gs- 1 ) 




0.5 


1.20(1.13-1.26) 


0.91(52) 


6.3X10" 11 


1.2X10 32 


3.0xl0 15 


3.1xl0~ 3 


1.0 


1.21(1.16-1.24) 


0.89(52) 


6.4X10- 11 


1.2X10 33 


2.9xl0 15 


1.6X10 -3 


5.0 


1.22(1.15-1.29) 


0.86(52) 


6.6X10" 11 


1.3xl0 33 


2.8xl0 15 


3.4xl0~ 4 



Table 7. The effect fixing fa to a range of values in the CWRK model during the spectral fitting. Details are given in S6.4 
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Figure 5. Results from spin phased-resolved spectroscopy of 
ASCA data. 0.4 < <p < 0.5 is the reference spin phase. The 
plots from the top are the spin folded light curve (included 
for clarity), the variation in the hard X-ray normalisation, 
the column density of the partial covering absorber and its 
covering fraction. 



normalisation, we go on to consider this explanation in more 
detail. 



6.1 Spin Pulse Modulation at Higher Energies 

It is evident from an inspection of the 8-25 keV spin-phased 
RXTE light curve (Figure 2) and the spin-phased normal- 
isation in the ASCA spectroscopy (Figure 5) that they are 
broadly similar. As this emission is expected to be at most 
weakly beamed, this indicates that the variation in inten- 
sity is caused by changes in visibility of the X-ray emitting 
region as the white dwarf rotates. The phase of maximum 
emission therefore corresponds to the phase of maximum 
visibility of the emitting region. There is possibly a slight 
phase shift between the two curves, with a clear peak at <j> 
= 0.2 in the RXTE light curve, and a more extended maxi- 
mum around 0.2 < <f> < 0.4 in the ASCA normalisation, but 
it is unclear given the uncertainties in the normalisations 
whether this is significant. Spin phases 0.2-0.4 also corre- 
spond to the phases at which the accretion region in the 
upper hemisphere is seen most close to face on (Potter et al 
1997, figure 9). This suggests that this is principle cause of 
the spin pulse modulation at higher energies. 



6.2 Spin Pulse Modulation at Lower Energies 

At lower energies, the maximum in the covering fraction of 
the absorber at <j> = 0.0 is the major effect on the soft X- 
ray light curve, while the local increase at <j> = 0.6 combined 
with the decrease i n no rmalisation at this phase to cause the 
second minimum (§4.2). Potter et al. (1997) find that accre- 
tion occurs preferentially along field lines which thread the 
disc ahead of the accreting pole (their figure 12). Material 
accreting along these field lines can therefore be identified 
as the source of the absorption at phase 0.0. Similarly, at (j> 
= 0.6, in the absence of absorption by the disc (cf § 3.4 ) , the 
local increase in covering fraction is likely to be caused by 
absorbing material along field lines which intersect the line 
of sight to the lower pole (Potter et al. 1997, figure 9). 

Figure ^ indicates that the change in absorption can be 
explained entirely by a change in covering fraction, with the 
column density of the absorber remaining constant. With 
reference to figure 9 of Potter et al. (1997), at cf> = 0.0, all of 
the field lines between the disc and the accreting pole will 
intersect the line of sight. At <j> — 0.2, those field lines to the 
leading part of the accretion region will no longer intersect 
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the line of sight, while by <f> = 0.4, only field lines to the 
trailing part of the accretion region will do so. If the column 
density through accreting field lines is ~ 5 x 10 22 cm 2 , the 
observed variation in covering fraction in Figure ^ can be 
reproduced. 

The accretion flow may be more finely fragmented along 
field lines feeding the leading edge of the accretion region. 
This would be expected because finely fragmented mate- 
rial is threaded by the magnetic field more easily than the 
larger denser inhomogeneities (e.g. Wickramasinghe 1988). 
From considerations of the packing fraction, the line of sight 
through more finely fragmented material is less likely to pass 
between gaps in the flow than in the case of the larger inho- 
mogeneities. This effect would reproduce the high covering 
fraction at cf> = 0.0, and its subsequent decrease towards <j> 
= 0.5. 



The unabsorbed s pec tral model from the analysis of the 
integrated spectrum (§4~lj) extrapolated for the energy range 
0.001-100.0keV gives the X-ray flux at 3.1 x 10~ 10 erg cm -2 
s _1 which, taking the distance to PQ Gem of 400pc (Patter- 
son 1994), gives a luminosity of 6.0xl0 33 erg s _1 and hence 
a mass transfer rate rate, M, of 1.4 x 10 16 g s _1 . This is typ- 
ical for IPs (Warner 1995). The resulting r A = 1-3 X 10 10 
cm or ~ 34Rwd may be too large given that an estimate 
of the distance to the first Lagrangian point is ~ 200 Kwd 
(Plavec & Kratochvil 1964) (the main uncertainty in rA is in 
the magnetic moment [iwd)- Nevertheless it does indicate 
that with a system inclination of 60° (Potter et al. 1997), the 
line of sight to the white dwarf surface is likely to be clear 
of the accretion disc at all spin phases, and this is unlikely 
to contribute to the cause of the covering fraction variation. 



6.3 Accretion Model 

The accretion scenario we are proposing fits neither the stan- 
dard occultation model of King & Shaviv (1984) nor the ac- 
cretion curtain model of Rosen, Mason & Cordova (1988). 
The orientation of the accretion region at spin phase max- 
imum is as predicted by the occultation model, but with 
absorption effects modifying the light curve significantly. In 
this it is similar to the "weak field/fast rotator" model of 
Norton et al. (1999) with the symmetry of the accretion cur- 
tain about the magnetic axis modified by the leading field 
lines preferentially stripping material from the inner mar- 
gin of the disc. The high (among IPs) magnetic field in PQ 
Gem and relatively high inclination and low dipole offset 
(~ 60° and ~ 30°, Potter et al. 1997) ensure that material 
on accreting field lines travels far enough out of the orbital 
plane to pass through the line of sight to the accretion re- 
gion, causing the observed absorption effects. In this it dif- 
fers substantially from EX Hya which possesses a magnetic 
field < IMG. 

6.4 Occultation by the Accretion Disc? 

Finally, we check whether the accretion disk can extend close 
enough to the white dwarf for it to have an affect on the X- 
ray light curves by obscuring the lower emission region. 

Using the Ghosh & Lamb formulation (Li, Wickramas- 
inghe & Rudiger 1996) the radius to the truncated inner 
edge of the accretion disc, ta, is given by 

r A = 0.52^(2GMr 1/7 M" 2/7 

where [iwd is the magnetic moment of the white dwarf. 

The magnetic moment can be estimated from the rela- 
tionship B = fi/r 3 . Using the fits to the RXTE data and 
the model of CWRK we found a best fit to the white dwarf 
radius, Rwd, of 3.8 x 10 s cm. Hence, with a magnetic field 
strength, B W d, of 15MG (Piirola et al. 1993, Vath et al. 
1996, Potter et al. 1997), [i W D = 9.6 x 10 32 G cm 3 . The 
accretion luminosity, 

L acc = GMM/Rwd 

where M is the mass of the white dwarf and L acc is emitted 
mostly in the X-ray energy band, enables estimation of the 
accretion rate, M. 



6.5 The size of the Accretion Region 

Using the radius of the white dwarf, the specific mass ac- 
cretion rate, m, (2.3 g cm _2 s _1 cf Table ^) and M de- 
termined above we can derive a fractional accretion area, 
/, of 9.0xl0 -3 . This is within the normal expectation of 
0.001;$ / <, 0.02 for an IP (Rosen 1992). However, from the 
RXTE data M =2.9 x 10 15 g s" 1 which implies a fractional 
area of only 7.1 xl0~ 4 . Refitting the RXTE data to the 
model in which m was fixed at 0.5, 1.0 and 5.0 g cm 2 did 
not give a significant adverse effect on the fit (x 2 =0.91, 0.89, 
0.86, respectively). Table |^ gives the implied fractional area 
and M for these m as well as the results from the spectral 
fitting. For a low specific mass accretion rate we find that 
the implied fractional area is consistent with the lower limit 
determined using previous observations. 

6.6 The Mass of the White Dwarf 

Previous determinations of the mass of the white dwarf 
in PQ Gem using an emission model fitted to Gmga data 
(Cropper et al. 1998, 1999) gave estimates > 1.1 Mq. In the 
case of the IP XY Ari, Ramsay et al. (1998) found that there 
was a good correspondence between the estimates given by 
this model and those from eclipse mapping. In our work it 
is found that the RXTE data gives estimates which are very 
much better constrained than those made with the ASCA 
SIS data. The estimate from our RXTE data of Mw d = 1.21 
(1.16-1.28) M & corresponds well to that of 1.21 (> l.O8)Af 
obtained with Ginga data. Although this appears to be un- 
usually high Ramsay (2000) found that the white dwarf in 
magnetic CVs were biased towards higher masses compared 
to isolated white dwarfs. 
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